Kondo screening cloud effects in mesoscopic devices 
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We study how finite size effects may appear when a quantum dot in the Kondo Coulomb blockade regime 
is embedded into a mesoscopic device with finite wires. These finite size effects appear when the size of the 
mesoscopic device containing the quantum dot is of the order of the size of Kondo cloud and affect all thermody- 
namic and transport properties of the Kondo quantum dot. We also generalize our results to the experimentally 
relevant case where the wires contain several transverse modes/channels. Our results are based on perturbation 
theory, Fermi liquid theory and slave boson mean field theory. 



I. INTRODUCTION 

In the last years, a large number of experiments have been 
realized in order to probe some theoretical predictions related 
to the Kondo effect. This has been made possible using small 
semiconductor quantum dotsi^ and also thin single wall car- 
bon nanotubes. 4 In these devices, it is expected that when the 
ground state of the quantum dot or of the thin carbon nan- 
otubes has a spin S = 1 /2 or equivalently when its last full 
level contains on average one electron and becomes degener- 
ate, the dot plays the role of the Kondo impurity when embed- 
ded into conducting leads. 

One of the major and still controversial issue in the Kondo 
physics concerns the Kondo screening cloud. The heuristic 
picture is that an electron in an extended wave function with 
a size of order £r- = Hvf /Tk (the Kondo length scale) sur- 
rounds the impurity, forming a singlet with it. The remaining 
electrons outside the cloud do not feel the impurity spin at 
low energies but rather a scattering potential, caused by the 
complex formed between the cloud and the impurity, result- 
ing in a 7r/2 phase shift at the Fermi energy (some recent ex- 
periments by Heiblum et a& have tried to measure this 7r/2 
phase shift but obtain results in disagreement with theoreti- 
cal predictions. 6 Nevertheless, the issue of the eventual rela- 
tion between the phase experimentally measured and the 7r/2 
phase shift is still controversial 7 ). The Kondo temperatures in 
the experiments on quantum dots are generally smaller than 
1°K and most importantly can be tuned via the gate voltage 
V g . Therefore, the Kondo length scale is expected to be very 
large. Typically with a Fermi velocity vp = 5.10 5 m/s, 
is of order 3 microns. Such low tunable Kondo temperature 
may therefore offer new opportunities to find evidence of this 
screening cloud. 

Recently, we have proposed that persistent currents in a ring 
containing a quantum dot may offer a way to probe d irectly 
this large Kondo length scaled See also Refs l9ll(ll Ill2ll3l 
Il4ll5ll 61). When the Kondo screening cloud becomes of or- 
der of the ring size or larger, we expect a crossover from a 
(large amplitude) saw tooth shape to a (small amplitude) sinu- 
soidal shape. Unfortunately, persistent currents experiments 
are generally delicate and very sensitive to disorder. In or- 



der to probe this fundamental length scale, it has been highly 
desirable to have some theoretical predictions based on con- 
ductance measurements which are more commonly used and 
better controlled. 

Possible devices to probe the Kondo length scale might be 
a quantum dot embedded into or connected to a finite size 
geometry. Several geometries where some finite size effects 
due to the large Kondo length scale may occur have been pro- 
posed recently: possible geometries include a quantum dot 
connected to a boxr^^ a mesoscopic ring 19 or a quantum dot 
embedded into short wires. 20 - 21 The device we want to study in 
detail in this paper consists of a quantum dot embedded in two 
short wiress (of length comparable to the Kondo length scale) 
which are then weakly coupled to some reservoirs. Such de- 
vice is schematically depicted in Figure 1 . One may wonder 
whether the conductance is sensitive to whether the Kondo 
cloud leaks into the reservoirs or is trapped inside the short 
leads. In the latter case, we may expect the conductance to be 
mainly dominated by the contacts between the reservoirs and 
the short leads. The answer is not so simple in the former case. 
Some partial results about this geometry have already been 
published elsewhere. 20 We have indeed studied a symmet- 
ric geometry, assuming the wires contain only one transverse 
channel, by means of perturbative calculations combined with 
a Fermi liquid approach when available. We have shown that 
when the bare Kondo temperature (defined by the Kondo 
temperature for infinite wires) is of the order or smaller than 
the energy level spacing in the wires, then the genuine Kondo 
temperature associated with the impurity can be very different 
from T 1 ^ and strongly depends on the local density of states 
seen by the dot being on or off resonance. In this paper, we 
would like to provide a more detailed analysis of this geom- 
etry and discuss the generality of the predicted finite size ef- 
fects. The plan of the paper is the following: in section 2, 
we present the tight binding model we want to study and give 
detailed derivations of these finite size effects using perturba- 
tion theory in the Kondo couplings. In section 3, we extend 
these results to wires which contain several transverse modes 
or channels. We note that the case of infinite length multi- 
channel wires has never been fully discussed to the best of our 
knowledge. In section 4, we show how these finite size effects 
affect quantitatively the transport properties through the quan- 
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FIG. 1 : Schematic representation of the device under consideration. 
en and e^ R control respectively the dot and wires gate voltage. 



turn dot. In section 5, we present complementary results valid 
at low temperature obtained with the Slave Boson Mean Field 
Theory (SBMFT). Finally in section 6, we discuss the main 
approximations made and the generality of our results. Three 
appendices with technical details are at the end of the paper. 



II. FINITE SIZE EFFECTS IN A BOX GEOMETRY 

A. Model Hamiltonian 

In this section, we want to study the following system: a 
quantum dot embedded in a quantum wire which is in turn 
connected to external leads by weak tunnel junctions. We as- 
sume that a gate voltage can be applied to the dot and also 
to the quantum wires. A related device has been proposed 
recently by Thimm et alM> where the quantum dot modeled 
by a standard Anderson Hamiltonian was coupled to all en- 
ergy levels of a finite size box with the same (energy inde- 
pendent) tunneling amplitude. The energy level spacing was 
assumed constant and of 0(1/ V) where V is the volume of 
the (3 dimensional) box. These two aspects differ consider- 
ably from the geometry studied here, since the electrons need 
to pass through the dot to contribute to the conductance. In 
fact, in the geometry studied in 1 17], the box is "side-coupled" 
to the dot whereas, here, the box embeds the dot. Moreover, 
the Non Crossing Approximation used in 1 17] might be ques- 
tionable in such geometry where several new energy scales 
em erg e compared to the usual Kondo model as was shown 
in 1 20]. We consider a one-dimensional tight-binding model 
which describes our device and has the Hamiltonian: 
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where the subscripts L, W and D stand for leads, wires and 
dot respectively. Here: 
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t l LW , with i = L,R (this notation will be used throughout the 
paper), are the hopping amplitudes between the different sec- 
tions of the geometry under consideration (see Figure^. The 
spin indices have been omitted in order to lighten notations. 

This Hamiltonian is a very simplified version of the geome- 
try of interest in this paper. We assume that the quantum wires 
are in the ballistic regime and do not contain any other impuri- 
ties (magnetic or not). We first consider the case of 1 channel 
where partial results have been reported elsewhere. 20 We will 
also discuss further a more realistic modeling of the quantum 
wires which typically contain 5 to 10 transverse channels^ at 
the end of this section. We ignore electron-electron interac- 
tions in the wires and leads, only keeping them in the dot. 
This point will be discussed in section lVTl 

We assume that the system is in the strong Coulomb block- 
ade regime, so that twD « —£d, U + eo, where eo < 0. 
Then we may eliminate the empty and doubly occupied states 
of the dot, so that Hp + %b = Hk gets replaced by a 
Kondo interaction plus a potential scattering term: 
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Here S is the spin operator for the quantum dot. 

The Kondo effect can be understood as resulting from a 
renormalization of the Kondo coupling constant, J, to large 
values at low temperatures. Perturbation theory is infrared 
divergent but the temperature acts as an infrared cut-off yield- 
ing a finite result which is accurate if the temperature is suf- 
ficiently high (T >> Tk). At low temperatures, a non- 
perturbative description is needed. This is provided by the 
local Fermi liquid description. 2 — If we imagine that J » t, 
then a spin singlet forms in the groundstate from the impu- 
rity and an electron in a symmetric orbital on sites ±1. The 
anti-symmetric orbital still remains available to conduct cur- 
rent so the system is roughly equivalent to the [7 = model 
and exhibits resonant conductance with the resonance tied to 
the Fermi surface. 9 

In the case of a closed ring, we showed earlier that the 
renormalization of the Kondo coupling is cut off, even at low 
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temperatures, by the ring circumference. In the present situa- 
tion, this renormalization would be cut off by the finite length, 
L = di, + dfi, of the quantum wires, if ti,w — 0. Essentially, 
if the Kondo cloud doesn't have enough space to form, then 
the growth of the Kondo coupling constant is cut off. The ef- 
fective Kondo coupling at the length scale L becomes of O(l) 
when L w £ K , Nevertheless, the situation is not so simple 
for small but finite t^w The growth of the Kondo coupling 
is no longer cut off by the finite size of the wire. Nonetheless, 
some noticeable finite size effects are expected to occur when 
L is reduced to a value of which corresponds to the 

screening cloud beginning to leak into the leads. 

B. Kondo temperature of the system 

Before computing any thermodynamic or transport proper- 
ties of the system, it is crucial to know the genuine Kondo tem- 
perature of our system, Tk, the main energy scale of the prob- 
lem. In order to calculate Tk, we first diagonalize the Hamil- 
tonian at J = 0, i.e. we diagonalize Hq = Hl + Hw + Hlw- 
Since the sites at j < and j > are completely decoupled 
in this limit, we can diagonalize the parts corresponding to the 
right and left leads separately. Let us first focus on the left part 



only. If t\ w — 0, then the wave-functions and eigenvalues of 
the left part (corresponding to — dz < j < — 1) are: 

ipL(j) = (1 fy/dl) sin k L . n j 
k L ,n = Trn/(d L + l); l<n<d L (7) 
e{kh.n) = —2t cosfc Lj „ + e w - fi. 

For non-zero t\ w , the spectrum of Hq (the left part of Hq 
is continuous). In order to study how the Kondo interaction 
renormalizes, we express c_i in terms of the eigenstates, c e 
of H^-hNl: 

c-x = / def L (e)c L , e . (8) 

We have normalized the operators cl £ so that {c^ e , cl i£ '} = 
S(e — e'). The left local density of states is defined by = 
|/l(c)| 2 and is normalized according to f_ 2t _ depL^) = 1. 
This local density of states can be easily computed exactly 
for our tight binding model. The calculation of pl(c) can be 
found in the appendix |A| The final result, assuming — 



PL{k) = 1 il sin2 fc ^ sinfc (9) 

Trt sin 2 fcz,(dL + 1) — cos fcsin(/cidL) sin(fci((ii + 1)) + 7^ sin 2 (fcic?i) 

where is related to k by 

— 2tcosk — (i = — 2tcoskL + (\ v — fi- (10) 

The fine structure and the properties of this local density of states has been studied in the appendix|X] For small t\ w , this local 

density of states, pl{z) = |/i(e)| 2 , has sharp peaks at the energies e„ = e\y — (J, — 2tcos[kL >n ], where the momenta „ are 
solutions of the equation: 



7 L sm(fc + fc L ) 7l sm(k - k L ) ., 

2k L (d L + 1) + arctan- 5 — - arctan- 5 — = 27m, (11) 

1 - 7£ cos(fc + fci) 1 - 72 cos(fc — k,L) 



with the constraint between k and k^ given by Eq. dlOi . For small jl, the solution of this equation is approximately 

k L ,n w nn/(d L + 1) + 0{~, 2 L /d L ). (12) 

Plugging this value of „ in (II Oi provides the value of the gate voltage necessary to reach such a resonance. The separation 
between two consecutive peaks is approximately w 2tt sin(fc^ „) /d^. The width of these peaks can be also evaluated as 

2(^ H/ ) 2 sin 2 (fc L , ra )sin(fc) 
t[d L + 1) 

We can deduce that the ratio of width to separation for a peak defined by fcz,, n is of order: 

x /a ~ l lw sin(fc £: »)sin(fc) 

OL,n &-L,n ~ —3 ■ (14) 

We will assume that this quantity is << 1, a condition already encountered for 7^ « 0.5. Obviously, similar results are obtained 
for the right part of Hq by changing L — > R. 
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The full Hamiltonian may be written in this basis as: 
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(15) 



where we have defined Jll — 2k| J, J RR = 2k r J, Jl R = 2klk r J. Similar definitions hold for the couplings V. It is worth 
noticing that these notations should not hide the fact that there is only one genuine Kondo coupling J and therefore only one 
Kondo temperature. Nevertheless, these notations are useful for transport properties as we be seen below. 

Using the Hamiltonian in the form of Eq. d!5l >. one can study how the Kondo couplings renormalize when the band-width is 
lowered from ±Dq (where Dq is O(i)] to ±D. At second order in the Kondo couplings, we obtain: 



Jij — > Jij + _ JikJjk 
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where k — L or R. The renormalization of the Kondo 
couplings is quite different depending on how far we lower 
the cut off, D. If D » A i n (i = L.R), the integrals in 
Eq. dl6l average over many peaks in the densities of states 
so their detailed structures become unimportant and we obtain 
the result for the usual Kondo model (see for example 1 25[). In 
particular, for a completely symmetric geometry, we recover 
the standard result J — > J[l + 2 Jp ln(D /D)} (J = J L l = 
Jrr = Jlr)- 

On the other hand, for smaller D, D << A;. n , the renor- 
malization of J in Eq. d 1 6t > becomes strongly dependent on 
the fine structure of the local densities of states and therefore 
of e\ w and e^- Let us analyze this in details in the next sub- 



sections. 



1. Fully symmetric geometry 

In order to simplify the picture, we first consider the com- 
pletely symmetric case. In this subsection, most of the sub- 
scripts L and R have been suppressed. Let us first assume that 
ew is tuned to a resonance of the density of states of width S n , 
Then the integral in Eq. dl6l gives a very small contribution as 
D is lowered from A„ down to 6 n so the renormalized Kondo 
coupling at the scale D, J e ff(D), practically stops renormal- 
izing over this energy range. Finally, when D < S n , the den- 
sity of states grows rapidly. To go further, we can approximate 
the local density of states (at the left or right part of the dot) 
by a sum of several Lorentzians of width S n : 
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example we have plotted in figure [2] the exact result and the 

approximate result for 7 = 0.5. We see that both curve agree 
2 




FIG. 2: The exact (plain blue curve) and the approximate (dashed 
red curve) local densities of states as a function of e. We have taken 
d = d L = d R = 49 and 7 = 0.5. 



remarkably well even for 7 = 0.5. They differ slightly in the 
height of the peaks and in the peak positions especially at the 
edges of the band. 



where k i>n = ■nn/{d l + 1) + 0(j 2 /di) (i = L = R). This 
constitutes a quite good approximation when ) 2 « 1. For 



We can express the result in terms of the change in J e f f (D) 
as D is lowered from A„ around the resonance peak at e = e n : 
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Je//(An) * ln( ^ )] ' 
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(18) 



where d = dj, = d^. When ew is tuned exactly 
on resonance (i.e. ew = £n)> the density of states ap- 
pearing in this renormalization is enhanced by a factor of 
(2isin 2 ki >n /[dsmk n S n ]) « t 2 /[t 2 LW sin 2 fc n ], fc„ being re- 
lated to fci.„ by J10I . This leads to a rapid growth of J e //(T). 

On the other hand, if ew is off -resonance then the density 
of states is small, of order ((5„/A 2 d) so the growth of the 
Kondo coupling is very slow at all energies D < A n . 

It is clear that this difference in the renormalization process 
between the two cases will strongly affect the Kondo temper- 
ature of the system. We define the Kondo temperature as the 
temperature where J e f / (T) becomes of 0(1). When J e f f (T) 
becomes large at T >> A„ then Tk is related to the bare 
Kondo coupling and bandwidth as in the usual case (with no 
weak links): T K » T K = D e~ 1/2Jpo . Furthermore, in this 
case, Tk does not depend strongly on ew- We may character- 
ize this case by T K >> A n or equivalently « L. The 
screening cloud fits inside the quantum wires and the weak 
links do not modify the Kondo effect significantly. 

On the other hand, suppose that T K << A n implying that 
J e //(A„) << 1. In this case Tk depends strongly on ew- 
If the system is tuned to a resonance then Tk will be slightly 
less than S„ : 
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On the other hand, if the system is off-resonance then Tk < < 
S n , the Kondo temperature reads 
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FIG. 3: The curves represent Tk = ,f(T K ) in a log-log scale keep- 
ing the same values for A n , S n for ew on resonance (plain curve 
which becomes almost flat at low T K <C 8 n ), and ew off resonance 
(dashed curve which drops sharply at low T K ). Both curves coincide 
atrS.>A n . 



Since the Kondo temperature, the main energy scale, enters in 
the calculation of all most all the properties of the system, this 
strong difference of behavior between the on or off resonance 
case will affect dramatically all transport or thermodynamic 
quantities. 



In this equation the subscript NR stands for non resonant. In 
general, for small values of tLw/t, we expect T^ R -C T§ 
(for experimental purposes, T^ R is almost 0). We have plot- 
ted the behavior of Tk vs. T k in Figure [5] for both ew on 
resonance and off resonance. This figure clearly illustrates 
the change of behavior when T K > A„ or T K -C A„. The 
curves coincides for T K ^> A„ as expected. WhenTj^- ~ A„, 
the Kondo temperature begins to be sensitive to the fine struc- 
ture in the local density of states which are signaled by the 
small oscillations of both T§ and T^ R (for the non reso- 
nant case, we have chosen the Fermi energy symmetrically 
between two resonant peaks) . Note that there is a regime 
(A/2 < T < A) where T^ R > T§. It corresponds to the sit- 
uation where the integrals in dl6l are dominated by two peaks 
for the off-resonance case and by one peak only for the on- 
resonance case. When T K -C A„ the off resonance Kondo 
temperature T^ R drops sharply at T K < A„ to very small 
values (<C S„). On the other hand T# also has a sharp drop 
at T K < A„ but then becomes almost flat and of order 5 n . 



2. Non-symmetric geometry 

If we consider now the general case of a non symmetric ge- 
ometry, 3 different situations may occur by tuning the two gate 
voltages e w and e^: (i) they are both on resonances, (ii) both 
off resonances, or (iii) just one is on resonance. In case (i) 
or (ii), the physics will be very similar to what we discuss in 
the fully symmetric case. A natural way to define the Kondo 
temperature is to evaluate D such that the ratios between the 
second order order corrections to the Kondo couplings and 

their bare values J^ll Jrr/Jrr, JlrJ Jlr ec l uals one - 
Since all the Kondo couplings are coupled together, we ex- 
pect this to occur at the same value of D = Tk, the energy 
scale of the problem. Provided the local densities of states are 
of the same order of magnitude, the discussion in the previ- 
ous section can be repeated. New features emerge for case 
(iii). Let us suppose that e w is tuned on resonance and e w 
off resonance. If T K is much larger than max(A L , A R ), we 
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obviously recover the bulk situation with Tk = T^- When 
Tj£ is smaller than the level spacing, then Tk will depend 
on the fine structure of the densities of states. The RG equa- 
tions d!6l > are then dominated by the resonance peak in the left 
lead. It means that the impurity is essentially screened in the 
left part of the device. Neglecting terms involving pu in the 
flow \\6\ . the Kondo temperature associated with this hybrid 
situation (we use the notation for this hybrid case) can be 
evaluated as before as 



with 



T, 



H 



' Do 



2(t£ w ) 2 sin 2 (fc„)/t 2 



(21) 



We therefore also expect the Kondo temperature to be of order 
6„ . Nevertheless, the strong asymmetry between the left and 
right part will affect transport properties. 



III. EXTENSION TO THE MULTI-CHANNEL CASE 

So far, we have modeled the short wires by a 1 dimen- 
sional tight binding model containing only 1 transverse chan- 
nel. Nevertheless, the semiconducting wires contain in gen- 
eral 5 to 10 transverse modes. In this section, we discuss ex- 
tensions of the previous results to take into account this multi- 
channel situation. Before considering the general case of finite 
length wires, we report for completeness on the more standard 
case of infinite wires connected to a quantum dot. 

A. Case of infinite wires 

1. Tight binding formulation 

This situation can be easily taken into account by introduc- 
ing in our tight binding Hamiltonian formulation a channel 
index (or band index) a — 1, . . . , N to the electron opera- 
tors which labels the electronic eigenmodes in the transverse 
direction. 

The tight binding Hamiltonian for infinite wires reads 
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Note that the spin indices have been omitted. It is straight for- 
ward to extend the analysis developed in the previous section. 
After a Schrieffer- Wolff transformation, the Kondo interaction 
reads: 
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In this section, we neglect the 
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direct scattering terms since they turn out to be unimportant 
for our discussion. In order to study the renormalization of 
this set of Kondo couplings, it may be more convenient to mo- 
mentarily suppress the distinction between left and right lead 
and to directly work with 2iV channels. This can be done by 



Cj,a+N, ~ * Ja+N,f3+N and J ( 



LR 

a,/3 



J, 



a+N,/3- 

We will switch between these two notations throughout this 
section. With these notations, the RG equations take a simple 
form 



Si 



Jay + Jaf J/3/3 / dt 



(25) 



where we have introduce pa (e) the local density of states seen 
by the dot in the channel 1 < (3 < 2N. Note that the equa- 
tion d25i follows from the fact the Kondo coupling J a p can 
be written as a product of a term related to channel a with a 
term related to channel (3. This product form is preserved un- 
der RG at least to lowest order. One may also introduce the 
dimensionless Kondo couplings constants 



A a ,/3 — \fPaPpJ a ,l3- 



(26) 



Since we consider infinite wires, it seems reasonable to 
assume that the densities of states are constant (p a = 
sm(kF,a)/^t m our tight binding model) in each channels 
such that the RG equations can be further simplified to 



dA 
dlogl 



= A 



0(A 3 



(27) 



where A = (\ a ,p) is a matrix of Kondo couplings. As usual 
the RG equations make sense provided < 1. All the 
couplings constant are driven to the strong coupling regime. 
Nevertheless, there is only one scale in the problem. This 
can be seen by noticing that the matrix A (which defines the 
bare values of the Kondo couplings) has only one non zero 
eigenvalue (= 7>A°). This is directly related to the fact that 
only one "effective" channel couples at the boundary to the 
impurity. Therefore, the Kondo temperature can be defined 
by: 



T K = Dcxp \-l/Tr{K 



,(o) 



(28) 



At T <C T^-, the impurity is screened by this "effective" chan- 
nel. It is straightforward to extend the analysis of Ng and 
LeeS (which relies on the Langreth description of the low 
temperature Anderson model 24 ) in order to prove that the con- 
ductance through a quantum dot symmetrically connected to 
two identical leads (containing several channels) reaches the 
unitary limit 2e 2 /h. In the more general case corresponding 
to two leads non identically connected to the dot, the T = 
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conductance in the Kondo regime reads 
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2e 2 
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with r£(e) = (i^n J 2 sin(A; F>a )/t. 



.a=l 



(29) 



2. Continuum limit formulation 

It seems also interesting to look at this multichannel case 
in the continuum limit. We have seen and concluded in the 
previous subsection that one "effective" channel couples at the 
boundary to the impurity. One may therefore wonder whether 
we can reduce the multi-channel problem to a single channel 
one. In order to formulate this model in the continuum limit, 
we may linearize the spectrum of the left/right leads around 
the 2N different Fermi points as: 



z L J R {x) =c a (±x) =e 



(30) 

V x > and where tp L ' represents a left moving electron 
in the left or right lead. Instead of working with left and right 
movers living on a semi infinite line we can work with left 
movers only on the infinite line using the fact that we have per- 
fectlyreflecting boundary conditions. The linearized Hamilto- 
nian reads: 



Hlin — 



N F 

1 E E IzWlJ^jtfxVl. 

= 1 i=L,R 



5 ^YftMiAx)+H.c.]\+H D , (31) 



with t % a — 2ismkp t at\ VD a , — 2isin(fc Fa ) and i labels 
the left and right side of the dot. Since the impurity couples 
only to one linear combination of operators ip l L a (0), it may 
be tempting to introduce a new field basis defined by a unitary 
transformation U l to the field ip l a defined by 



Vi 



(32) 



The other set of operators rf a with a > 1 can be built in 
such a way for the matrix U = (u) a .j3 to be unitary by the 
Schmidt orthonormalization process. Nevertheless, this uni- 
tary transformation leaves the kinetic part diagonal provided 
= v F Va. If this condition is met, the Anderson impurity 
model reduces in this new basis only to one effective channel. 
A similar conclusion was recently reached independently by 
Cho et al. 22 considering several independent leads (containing 
one transverse channel). Since at low temperature, transport 
depends only on the Fermi level properties, the conductance 
through the dot will be dominated by one effective channel. In 



particular, for a symmetric geometry, we recover that the con- 
ductance may reach at most the value 2e 2 /h corresponding to 
the unitary limit. Nevertheless, for several transverse modes 
in a wire the condition v F = vp is most likely not satisfied 
and the introduction of the unitary transformation U may be 
not very helpful since it does not leave the kinetic energy di- 
agonal. (One can still argue that the difference between Fermi 
velocities is usually irrelevant in ID systems and the conclu- 
sions obtained above should hold.) Therefore, when the Fermi 
velocities are different, we need to perform some RG analysis 
analogous to what was developed in the previous section. In 
this formulation, the Kondo Hamiltonian reads: 



2JV 



H K = 



a, 13=1 



(33) 



where \ a< p has been defined in J26l >. Note that we have used 
the compact notation with 2N channels in order to avoid the 
distinction between left and right wires. Starting from Eq. 
J33l >. one can obtain in a similar manner the RG equations 
given in Eq. J27l > 



B. Case of finite size wires 

In this section, we extend our previous results to take into 
account the fact that the wires are now finite and contain sev- 
eral channels. A possible Hamiltonian which extends Eq. Q 
and describes this situation reads 



H — Hl + H\\/ + H 



H 



LW 



H 



WD 



He 



(34) 



with 



H L = -t 



— (II —2 00 

E + E 

j— — 00, a j=d_fj+l,Q 



E +E 

3=—dt,a 3=1, a 



( C 3,<x C 3+h<* + M-3 5 ) 



+4 



E 



E 



(36) 



Hlw 
Hwd 



j=-d L ,a j = l,a 



a,0 



-*£W,a/3 c L,a C <*fl+l,/3 + H - c -) (37) 

- E(*W c -i,a c o + twncla,* + ff.c.)(38) 



(39) 



Hp> is the dot Hamiltonian and remains unchanged compared 
to \l\ Note that the tunneling amplitudes between the wires 
and the leads mix together the channels. Allowing channel 
mixing will not change much our discussion as we will see 
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further. The different bands have dispersion relations 

£ a (k) = — 2tcosk — n a 



-2icos Ul ■ 
-2t cos kn 



Mo 



(40) 



1. Without channel mixing 

Before treating the general case of a matrix of tunneling 
amplitudes between the wires and the leads, we first consider 
in this subsection the simpler case where the matrices t^L are 
diagonal in the channel index. Such an approximation is valid 
when the wires and the leads are geometrically similar (i.e. 
with the same number of tranverse modes) and the tunneling 
barrier separating them has no dependence on the transverse 
direction. 

The fact that thw.u ^ t does not modify much the dis- 
cussion made in section ITlI Al A Schrieffer- Wolff transfor- 
mation leads to the same Kondo Hamiltonian given by Eq. 
{23). The main difference is that the local density of states 
seen by the quantum dot acquires some non trivial energy de- 
pendence. Therefore, the RG equations corresponding to Eq. 
d25i are still valid. The density p a (e) generalizes Eq. © with 
7l — * 7a = t W L,a/t> k L — + k a , fi — > n a and corresponds to 
the LDOS in channel 1 < a < 2N (within this compact nota- 



tion t$r Ll 



twL,-y, t WLn = t WLn+N y~/ e [i,iv])- Note 

that there is still only one scale in the problem, the Kondo tem- 
perature Tk = D cxp[— /7>(A(0)], which depends on the lo- 
cal density of states seen by the impurity, more exactly on the 
total tunneling amplitude T(e) = irJ2 a ^wd a P«( e ) ( see Eq. 

The properties of p a are therefore similar to those of pl(c) 
studied in section 111 Bl Let us for example analyze the total 
left tunneling amplitude seen by the quantum dot defined by 

N 

t\vD.aPa( e )- Since the \i a are different, 



a=N+l 



J 



In this basis, the Kondo Hamiltonian reads: 



the position of the peaks in the LDOS pL, a are in general 
different. It implies that the average level spacing in the 
quantity Ti,(e) is now of order Aj^ = ttvf/NcIl ~ Al/N. 
Moreover the peaks have now different widths depending on 
t WL , p a and e w . The analysis developed in section HlBl 
extends straightforwardly: when T^- becomes of order of or 
smaller than Aj^, finite size effects appear in a completely 
similar way. It also implies from an experimental point 
of view that it is preferable to work with thin wires with 
a few transverse channels in order to reach the condition 
Tjf < A/N. It is also worth mentioning that while in the 
1 channel case, speaking about the ratios between A/T^ 
or £k/cI were equivalent, it is not the case here. We can no 
longer simply compare the Kondo screening cloud size to 
the length of the wire. Nevertheless, we can still define an 
effective length d e ft — Nd which can be compared with 

2. With channel mixing 

Suppose we now work with a more general model than \35\ 
with a tunneling matrix t^ L a „ between the left lead and the 
left wire (and the same in the right lead). We may express 
C-i. a in terms of the eigenstates, C£ :7 (e) of the left part of 
Hl + Hyv + Hlvv + Hq'- 



/2i-ju ^ 
de^/L,a,7(e)cL, 7 (e)- (41) 
-2t-jU 



The local density of states in the channel a at site —1 can 
be then defined by p Q (e) = |/i, Q , 7 (e)| 2 and is normal- 
ized accordingly. This LDOS has a level spacing of order 
A ~ Ttv^jdh provided all the elements of the tunneling ma- 
trix are small compared to t. It implies that the resonance peak 
positions of |/l jQ . 7 | are almost similar V 7 given by the ones 
of the isolated wire (up to a small shift of the order of the peak 
width which depends on the tunneling matrix). 



H 



K 



a,6e 7,6 
tLRJ 



(42) 



where the couplings have been defined in Eq. (I24> . One 
may define new tunneling amplitudes vt = ^ Q fL.a,-yt l WD a 
(with i — L,R) and a new set of Kondo couplings: 



XUe,e') = {viy{e)vl{e')/i D . 



(43) 



In this new basis, the problem takes a more familiar form 
with a 2N x 2N matrix of Kondo couplings (if we again 
switch to the compact notation where no dictinction between 



left and right is made): 



H K = 



2N 



dede' Ksc\(e)^ ■ Sc 5 (e') (44) 
7,5=1 



The RG equations have already been derived and read in this 
notation: 



A76 ^ ^7(5 + -^75 J]] J de 



(45) 
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The coupling A„„ (e) has some non trivial energy dependence, 
since its bare value \ vv {e) oc Y^ a \fa,v(e)twD, a \ 2 - The 
average resonant peak spacing of ^ Q \f a ,t>\ is 0(Tivf /Nd) 
and the position of the peaks are almost independent of v. 
Therefore, in order to calculate the Kondo temperature of 
the system, we need to compare to the energy scale 
0(hvp/Nd) = A/N as in the non channel mixing case. To 
summarize, for weak tunneling amplitudes between wires and 
leads, the mixing of the channels just affects the fine structure 
of the local density of states seen by the quantum dot (i.e. the 
width of the peaks) but not the resonant peak positions (and 
therefore the level spacing) which is almost determined by the 
diagonalization of the isolated finite size wire. We note that 
this conclusion is quite general and should also be valid for 
the more realistic experimental situation of finite size wires 



connected to 3D reservoirs. 



IV. STUDY OF TRANSPORT PROPERTIES 

We apply in this section the results obtained in the previous 
section in order to consider transport properties, in particular, 
the conductance through the device under consideration. As 
for thermodynamic properties, a perturbative calculation is in- 
structive. We expect it to be valid at sufficiently high T ^> Tk 
when the renormalized couplings are sufficiently small. The 
linear conductance (for t l LW ^ 0) at cubic order in J, V is 
given by: 



G(T) 



e 2 f 3 

— 7r 2 / dep L (e)p R (e)[-dn F /de]-J LR [l + J LL I L (e) + J RR I R (e) 



e 2 f 

^47rV L 2 fl / dep L (e)p R (e)[-dn F /de) 



(46) 



r 



with 



PL(e') 



(e'-e) 



and a similar expression for I R . PP stands for principal part 
and riF(e) is the Fermi distribution function at temperature T. 
These auxiliary functions leads to the usual logarithmic cor- 
rections. The derivation of this formula is sketched in the 
appendix 151 We have not written the other contributions at 
cubic order since they do not involve logarithm divergences 
(and equal zero for particle-hole symmetry) and therefore do 
not renormalize in the infrared limit (see appendixlBl. 

Notice too that the integral in the expression of Li (e) (i = 
L, R) depends on the local density of states pi(e). 

Let us focus on the second order terms in Jl R and Vlr, and 
ignore, for the moment, the corrections of higher order. We 
must distinguish 3 regimes of temperature resulting simply 
from the fact that (—drip/de) has a peak with width of 0(T). 
If T > > A l n , then the integral in Eq. J46> averages over many 
peaks of Pi(e) so that G is approximately independent of e^: 



where po = sin ftp /wt is the average local density of states. 

When 5 l n « T « A^, the conductance depends 
(1 — 2np(e')) (47) strongly on the gate voltages e z w . We then encounter three 

possibilities: i) both e l w are tuned to a resonance peak, (ii) 
both are far from a resonance peak, (iii) only one of them, let 
us choose e^pr is tuned to a resonance peak. Let us detail these 
three cases: 

(i) If both e]y are tuned to a resonance peak such that 
e^r — p + 2tcos[fci,„] and e{y — p + 2tcos[fci iin ], then 
the integral in Eq. J46l > is dominated by the peaks located at 
kh.n ~ Kn/{d,L + 1) and k R , m w irm/(d R + 1) (notice that 
the integers n and m are in general different). In this regime 
of temperature, we can approximate (—dnp/de) « 1/4T. 
Moreover, we can neglect the other peaks except the resonant 
ones in the LDOS: 



2sin 2 (/c Li „) 



G 



— ^p G f[?,Jl R /A + Wl R ] (48) 
7rn, 



and we have a similar expression for irp R (e). Performing the 
integration we obtain for the on resonance conductance (sub- 
script R): 



C R m - e " ( 3J £fl/ 4 + W Ir) ^sin 2 (fc L:K ) sin 2 (fc fl , m ) 
[ ' ~ ttTi Td L d R # + 

- —(3J 2 /A + W 2 ) 7rtsin 2 (fc LiTl )sin 2 (fc rt . ; „„ 

~ nh { LRl LR> 2Tsink [(t^d^m^k^) + (t* w ) 2 d L S m 2 (k R>m )] 



■R,m) (49) 
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For a completely symmetric geometry, this expression reads: 

G R (T)^ __-M^__A2-^ M1 r (/,,,, I ,50) 



7r?i 2Tdt 2 LW sin A: 



with fe i)n = fc Lj „ = fc_R jm . 

(ii) On the other hand, if both e l w are far from a resonance peak (compared to T) then the LDOS is almost constant and may be 
approximated as 7rpi(e) w 2sm j^ 1 "^ A similar expression for /9r(c) can be written. Therefore the non resonant (subscript 
NR) conductance reads: 



iNRfrw* „ e2 tt r2 ,i, iV s , sinfc £;K sinfc^ m gjg 



G ^«(T) » ^(3Jl fl /4 + 4F L 2 fl ) ( ^ )2(t ^ )2 ^e'" Sln2 Sln2 < 51 > 
We also give the expression for the completely symmetric geometry 



e 



2 i a i /it/ 2 sm (^i,«) sm ^ 



G ™(T) « -(3Jl R /4 + 4 ^^™riJ:. (52) 
We immediately notice that G NR (T)/G R (T) = 0(t 6 LW /Tdt 5 ) < 1. 

(iii) Finally, if only is tuned on resonance, then the conductance associated to this hybrid (H) situation reads: 

e 2 , o o , sin /cr „ sin ko „5 R 

* ^(3^ + 4^ )z | ^ • (53) 

The conductance G H in this hybrid case is also small compare d to G fl : G H {T)/G R {T) = Q(tj w /Tdt 3 ) < 1. 



In the ultra-low temperature regime, T << S l n , one may 
evaluate the conductance in a similar way. by approximating 
npi(e) by 2siv?(ki^ n )/di8' l n . When and are both on 
resonance, we can approximate irpi{e) by 2 sin 2 ^^)/^^. 
The conductance then reads 



G{T) 



(3J 



Win 



)t 2 



ttTi 



\ i LW i LW 



) 2 sin 2 k 



(54) 



L and 



The conductance is still given by Eq. \5 1\ when e w anu c w 
are both tuned off resonance for T < S n . Finally, when only 
e w is tuned on resonance, the conductance reads: 



G(T)^—((3J 2 LR /A + Wl R ) 

TTIL 



sin 2 k Rtn 



(t R w) 2 



(55) 



We have used perturbation theory in order to evaluate the 
finite temperature conductance in all different situations. We 
have seen that the on-resonance conductance (where both 
wires are tuned on resonance) is far larger that the conduc- 
tance in the other situations. This conclusion is valid only 
where perturbation theory applies. These approximate formu- 
las certainly break down when they do not give GttTi /e 2 < < 
1, due to higher order corrections in J and V. So our approx- 
imate formulas will certainly break down before T is lowered 
to 8 % n unless J « {t l LW ) 2 ft, a condition which might typi- 
cally not be satisfied. When these formulas apply, we clearly 



see that the conductance is much larger when e w are tuned on 
resonance. 

However there is another, more interesting reason why 
these formulas can break down at low T, namely Kondo 
physics. The cubic correction in Eq. J46i contains a In T term 
which essentially replaces Jl R by its renormalized value at 
temperature T, J e J R (T) following d!6l >. We expect that this 
will remain true at higher orders. 

Now consider the behavior of the conductance as a function 
of T and e' l w in the two cases. In the case << d,-, we may 
calculate the conductance perturbatively in J e ^ R (T) at T >> 
T# and using local Fermi liquid theory for T « 1%M For 
T >> Tft, we obtain Eq. J48> . essentially independent of e l w . 
On the other hand, for T « T^, the conductance reduces 
to that of two reservoirs connecting each other by a contact 
whose conductance is 



G 



w 



2e 2 /h {2t^ L t R L ?K{twL? 



(4l) 2 ) 2 



(the low temperature conductance of the quantum dot in the 
Kondo regime). It corresponds to our original model with 
U = 0, € c i = 0, and some effective length di ~ di. (di 
can be somewhat reduced from di by an amount of order 
A general but less readable formula for the transmission asso- 
ciated to such non interacting geometry is given in appendix 



When t& L 



*wl and e w 



ew, the quantum dot 



reaches the unitary limit and the conductance reduces to that 



11 



of an ideal wire of length L = <Il + cLr. The expression of such conductance reads 
I 



G(T) = P — 
no, 



2t sin k dk 



4:(tt LW ) 4 sin 2 k d sin 2 k 



A 2 + B 2 



with kd defined by cos kd — cos k — ew/2t, and 



~de~ 



■(*)] 



A = 
B = 



sin[fc d (2L + 2)] - (7! + j%) cos k sm[k d (2L + 1)] + 7^7! cos 2k sm[2k d L] 
-(jI + 2jj,)smksm[k d (2L + 1)] + 7 2 7^ sin 2k sm[2k d ]L 



(56) 



As T is lowered below A l n the conductance for a symmetri- 
cally connected quantum dot develops peaks with spacing of 
order A' l n /2 (i = L, R). This is the spacing of peaks in the 
density of states of a wire of length L = d^ + d^, contain- 
ing no quantum dot. It is half the spacing in the density of 
states of the model with J = 0, discussed above. Initially, as 
T is lowered below A„, the peak width is of 0(T) and the 
peak height is of 0(2e 2 A n t 2 VL /hTt 2 ). As T is lowered be- 
low S n the peak width becomes of 0(5 n ) and the peak height 
becomes of 0(2e 2 /h). For a non symmetrically connected 
quantum dot (t WL 7^ t^- L ), tne halving of the period occurs 
when the Kondo quantum dot has a large enough conductance. 
In particular for two symmetric wires, we can prove using Eq. 
iA6\ that the number of resonant peaks doubles when the con- 
dition 



\Rd\ < 



2\R L 



w 



1 + \Rlw? 



(57) 



LWJ 



is satisfied. In this equation, we have defined Re> = (t 
&) 2 /(('w) 2 + {tfw) 2 )i *e reflection probability of the 
dot at low temperature and Rlw the reflection probability of 
the left or right weak link. 

On the other hand, when £k » di,, dn, the dependence 
of conductance on T and ew is very different. As T is low- 
ered below A l n the on-resonance conductance starts to grow 
both because of the single-electron effects reflected in Eqs. 
( I49> and J54i and, eventually, when T < S n because of the 
growth of J e ff(T). However, off resonance the conductance 
stays small, given by Eq. j5U at least down to temperatures, 
T « S n of 0(Tg R ), given by Eq. In the tempera- 

ture regime T^ R << T « A„, the conductance has peaks 
with spacing A„ reflecting the fact that J e f f (T) is small, off 
resonance. In Figure@]we sketch the conductance versus ew 
in the two cases T$- > A„ (dashed style) and T£ < A„ 
(plain style) reflecting the halving of the period between the 
two curves. Note that we have used perturbation theory to plot 
the plain curve explaining why the amplitude is small. Nev- 
ertheless, at lower temperature, it is more difficult to calculate 
the on-resonance conductance both because of the breakdown 
of the perturbative result of Eq. ( I49> . J54t due to single elec- 
tron effects and because it appears considerably more difficult 
to extract unambiguous predictions from local Fermi liquid 
theory. Nonetheless, as we will check in the next section us- 




e w /A 



FIG. 4: Conductance in a symmetric geometry as a function of ew at 
fixed A n , T and S n for both cases £a' 3> g(l = da (plain style) and 
£k <C dL (dashed style). We have chosen T ~ A„/10 3> 5 n ~ 
A n /100. The curves have been represented on phase but in general 
they are expected to be shifted, the shift being difficult to determine. 



ing slave boson mean field theory, it is very reasonable to ex- 
pect a conductance of O(l) on resonance at T < S n where 
J e ff(T) is O(l) on resonance. Off resonance we can show 
rigorously that the conductance remains small since J e / / (T) 
remains small there and so do the single electron corrections 
to Eq. J5 II . The behavior of the conductance at very low 
T < Tft R in the case £k >> d will be also analyzed using 
the slave boson mean field theory. 

In Figure [5] we have drawn schematically the conductance 
on resonance as a function of temperature for two different 
bare Kondo temperatures T^- 1 3> A n and Tjj 2 &n ^ 
A n , using the perturbative formula given by Eq. d46b and 
the Fermi liquid picture valid for the first case only. For the 
first case, the conductance has a plateau which corresponds 
to the quantum dot being screened and the e integral in Eq. 
j46t averaging over many peaks. The conductance reaches 
2e 2 /h only when T <C S n , Conversely, in the second case, 
the conductance remains small til T w S where the Kondo 
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FIG. 5: Conductance as a function of temperature in a symmetric 
device (assuming both e l w is on resonance) for both cases A n <C 
T%,i (right blue curve) and A„ > 8 n > t£_ 2 (left red curve). The 
curves in plain style correspond to the perturbative calculations plus 
the Fermi liquid result for the first case only. We have schematically 
interpolated these curves (dotted lines) where neither the perturbative 
nor the Fermi liquid theory applies. 

coupling becomes strongly renormalized (see Eq. dl8» . We 
may expect a very abrupt increase of the conductance in this 
regime as schematically depicted in Fig [5] Notice that for this 
choice of 2 , the renormalized Kondo temperature T^ 2 is 
actually enhanced and of order 5 n . These different behaviors 
lead to different shapes of the curves. 

V. SLAVE BOSON MEAN FIELD APPROXIMATION 

We will use in this section the Slave Boson Mean Field The- 
ory (SBMFT) in order to confirm qualitatively our analysis 
in the previous section. 26 This approximation seems to well 
describe qualitatively the behavior of the Kondo impurity at 
low temperature T <C Tr- when the impurity is screened. It 
has recently been used by Hu et al. U) to calculate the persis- 
tent currents in a metallic ring containing a quantum dot and 
their results nicely confirm our perturbative and renormaliza- 
tion group calculations. This method assumes U = oo where 
the impurity operator can be written as dj. = /J b, where the 
fermionic operator f a and the bosonic operator b describe the 



singly occupied electron and holes states respectively. Fur- 
thermore, the constraint ¥ b + ^2 flf a — 1 has to be im- 

a 

posed. In the mean field approximation, the boson operator is 
replaced by a c-number bo. The constraint is implemented by 
using the Lagrange parameter Ao- Therefore, in the mean field 
approximation, the Hamiltonian is the same as in Q with the 
following changes 

H WD = -boitfyndj + t^jjfld + H.c.) (58) 
Hd = e / t / + A (6 2 1 -l), 

where we have defined eo = £_d + Ao. 

Within this approximation, the Hamiltonian reduces to a non 
interacting system. The values of Ao and bo can be determined 
self-consistently by minimizing the free energy of the system 
defined by F M f = — jjlogZ + A (&q _ 1) with 

Z = Y[(l + e-^"- E ^) 2 (59) 
fc 

with Ep the Fermi energy and ej. the energy eigenvalues of 
the mean field Hamiltonian. Once bo and Ao are found, we 
may directly apply the Landauer formula to find the conduc- 
tance at finite temperature of the whole system. This is made 
possible because the mean field Hamiltonian is a non interact- 
ing one. We also want to emphasize that it is crucial to in- 
corporate from the beginning the reservoirs when calculating 
the numerical parameters (here bo and Ao) characterizing the 
quantum dot, especially in the limit dp, <1r <C Note that 
in 11711 . the quantum dot retarded Green function was calcu- 
lated numerically using the non-crossing approximation with- 
out incorporating the reservoirs (i.e. by taking t l WL = 0). 
This might be a crude approximation especially in the limit 

The mean field free energy is conveniently expressed as: 

DO 

Fmf = ~ J den F (e)Im(\nG r f (e)) + \ Q (b 2 o-l), (60) 

-Do 

where np{e) = 1/[1 + cxp(/3(e — fi))] and G r j is the retarded 
Green function at the impurity. Using the approximate expres- 
sion dl7> for the left and right local density of states, the free 
energy reads: 



DO 



Fmf = — 

IT 



de n f (e) arctan 



-Do 



( E 6i(*U) 2 ^i(e) 

i=L,R 

t-to-Eb 2 o(t l w D ) 2 ZUe) 

i 

i 



Hb 2 o - I), 



(61) 



where 



2 dl 

n— 1 v 



) 2 + s 2 

i) ' u i,i 



(62) 



n=l 



sin 2 ki 



(e - £»,n) 
(e - e*,„) 2 + S. 



2 • 

i . n 



(63) 



13 




CM 

O 



0.4 - 



0.2 - 



-rO,NR \ \Y 
'k,2 A V 


'K,2 


\ 

\ 

\ 

\ 

\ 
\ 
\ 
\ 

\ 

-y- 


\ 

H - 


'K2 i 

i . , ' i i 


* 



1°0" 7 10"° 10" 3 10"" 10"" 10"' 10" 

8 a 



T 



FIG. 6: Conductance from SBMFT as a function of temperature for 
two bare Kondo temperatures T^.i ~ 3.6A and Tr',2 ~ A/10. 
For the latter case, we have considered the 3 cases depending on e' w 
being both on resonance (corresponding to T% '%), both off resonance 

(corresponding to T^' 2 ) or the hybrid situation (corresponding to 
T%"). Here 7 = 0.3. We have not plotted the off resonance case 
because T^ H seems very small in this case and our SBMFT has 
convergence difficulties. 

By minimizing Fmf with respect to Ao and bo, we obtain a 
set of two self consistent equations which can be solved nu- 
merically by iteration. Since the formulation of the problem 
is non-interacting, the conductance can be expressed using the 
Laudauer formula 

where T(e) is the total transmission probability through the 
non interactive device in the SBMFT, the quantum dot is mod- 
eled by a simple resonant level with U — 0, t l WL — > bot % WL 
and eo = + Ao. The expression for T is given in the ap- 
pendix B by the Eq. (ICli . 

A. Analysis of the symmetric geometry 

As in the previous section, we have fixed d (therefore the 
level spacing A) and considered different bare Kondo temper- 
atures corresponding to both limits T^- <C A and T^- ^> A, 
where A is the level spacing at the Fermi energy inside the 
wires. In the following, we assume that t^ D — t^ D , 

Ilw = 1 lw = 7* and d L = d,R = d. 

We have plotted the conductance as a function of tempera- 
ture for two different bare Kondo temperatures T% 1 w 3.6A 
and T° 2 » A/10 for 7 = 0.3 (Figure^ and 7 = 0.5 (Figure 
0. For 2 , we have considered three cases: both left and 
right wires are on resonance, both wires are off resonance and 
finally only one of them is on resonance. Such a distinction 



FIG. 7: Same as Fig. |6|but with 7 = 0.5. Note that here the 3 cases 
are plotted for T£ j2 w A/10. 



is not worth making for Tjj 1 since all curves are very similar 
since the conductance does not depend of the fine peak struc- 
ture of the LDOS. As previously emphasized, the SBMFT is 
a low temperature method complementary to our perturbation 
calculation. Therefore, we have stopped the curves at T w Tk 
(not the bare one !) where the method failed or is meaningless. 
Moreover, some convergence problem may occur in this limit 
particularly for the case -C A where the Kondo couplings 
get very strongly renormalized on a very small energy scale. 
Another more powerful numerical method like the Numerical 
Renormalization Group (NRG) carried out in Ref. 12111 is re- 
quired to estimate the conductance for the full crossover from 
low to high temperature. 

Let us analyze the results now. In figure |6j the curve 1 
looks like the one predicted in Figure [5] It has a plateau 
around T = A separating the low temperature where the 
unitary limit is reached from the high temperature perturba- 
tive regime (almost not shown here). The curves for T^J 2 is 
different and does not exhibit such an intermediate plateau. 
On the other hand, when the temperature is increased, it sud- 
denly shoots down abruptly at T w T§ 2 w 0.002 ps 3(5. 
We quickly go from the low temperature regime to the high 
temperature regime where perturbation theory applies. Note 
that this sharp drop is also associated with the limit of conver- 
gence of the SBMFT in this situation (because the variational 
parameter bo becomes very small). This sharp change of be- 
havior in the conductance was predicted by the perturbative 
analysis (see Figure |5jl and is due to the strong renormaliza- 
tion of the Kondo coupling when the temperature reaches the 
order of the peak width. We have not plotted the conductance 
for the off resonance case where our result indicates a Kondo 
temperature T^f * 10 -7 with an abrupt change too. The on 
resonance- off resonance curve exhibits a similar behavior as 
the on resonance one but its value at low temperature is very 
small due to the asymmetry. For 7 = 0.5, we still observe an 
abrupt decrease of the conductance for 2 in the three pos- 
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FIG. 8: Conductance as a function of temperature in a completely 
asymmetric geometry with only one short wire. We took jl = 0.5 
and 7n = 1. We have considered two bare Kondo temperatures 
T£ %1 « 3.6A and I% <2 ~ A/10. 



FIG. 9: Same as above. We have considered three bare Kondo tem- 
peratures (see the labels) which all correspond to <C A. Note 
that in this figure A w 0.1 and 5 ft! 0.006. 



sible situations. On the other hand, the difference of behavior 
between T® 1 and T^'^ has almost disappeared . In particular, 
the plateau around T s» A for 1 is now replaced by a shoul- 
der. The only way to identify the situation <C A relies on 
the dependence on the local densities of states being on or off 
resonance. For this purpose, we note that the behavior for the 
off resonance curve defined by the label Tj^ R is very similar 

to the on resonance one defined by Tg 2 , but translated toward 
lower temperature by roughly two order of magnitudes. 

From this difference, we can immediately infer that the con- 
ductance function of the gate voltages = e™- exhibits a 
doubling of the period (compared to the case T£- 3> A) in 
the window T^f < T < 1% v Instead of varying two 
gates voltage, we may imagine that we fix on resonance 
and study the conductance as a function of only (which is 
more convenient to realize experimentally). Suppose the tem- 
perature is also fixed at T ss 5 (see Figure^). The curve corre- 
sponding to Tg- ^> A will exhibit resonance peaks of height 
« 2e 2 /h with the level spacing A/2 as depicted in Figure 
0] On the other hand, the curve corresponding to <C A 
will exhibit a succession of large resonant peaks(of height 
« 2e 2 /h) and small resonant peaks (of height w 0Ae 2 /h), 
clearly different from the case ^ A. 



B. Analysis of the complete asymmetric geometry 

One may wonder what happens if we connect only one 
short wire, the left one for example, to the quantum wire, 
an a priori easier situation to realize experimentally. In our 
model, it corresponds to = 1, In the figure |8j we have 
therefore plotted the conductance as a function of tempera- 
ture for the two Kondo temperature 1 = 3.6A « 0.36 and 



Tjr 2 = A/10 « 0.01 (defined for infinite wires). Except for 
the shoulder around T as A for the curve Tjj x , the curves are 
very much alike. For the geometry with two wires, the best 
signature of the behavior Tjj -C 8 was the difference between 
the on- or off- resonance situations. Here this difference is 
not very dramatic (except for a sharper behavior for T*b 2 ). 
The effective Kondo temperatures for both cases are around 
Tk 2 — 0.001 (slightly less for the off resonance case). It 
implies that the renormalization of the Kondo coupling is es- 
sentially dominated by the bulk part which does not contain 
any wire. In other words, the screening cloud builds mostly in 
the right lead. Nevertheless, it may not be always be so. In- 
deed, we have plotted in Figure [9] the on-resonance conduc- 
tance for three different Kondo temperatures TE 1 = 0.1 A, 
t k,2 = 0.005A and finally T°- 3 = 3.2 10~ 4 A. The first 
Kondo temperature corresponds to a situation studied in Fig- 
ure [8] where the screening is mainly dominated by the right 
lead. In other words, for such bare Kondo temperature, we 
do not expect much variation of the conductance with eh,. 
The other two bare Kondo temperatures have been chosen 
such that T^ 2 ^ 3 <C 5. For such choice, we expect that 
the renormalization of the Kondo couplings will be domi- 
nated by the on resonance peak in the left LDOS. This cor- 
responds, as mentioned in the previous section and exempli- 
fied here, to an enhancement of the Kondo temperature. From 
the conductance curves, one may infer a cross-over tempera- 
ture which can interpreted as an effective Kondo temperature. 
One could approximately evaluate Tk,2 ~ 0.002 > 2 and 
Tk.2 ~ 0.0005 3> T% 3 . In such cases, one expect a strong 
dependence of the conductance with e^. Indeed, when is 
off resonance, the renormalization of the Kondo coupling will 
be again dominated by the right lead with an effective Kondo 
temperature « Tjj 3 . It implies that the conductance 

as a function of should exhibit peaks of spacing A l for 



the temperature range T^f <C T <C T^ 3 . To summarize 
the analysis of this completely asymmetric geometry, finite 
size effects associate with the formation of the Kondo screen- 
ing cloud may occur only when the bare Kondo temperature 
is less that the resonance width. It implied very low Kondo 
temperatures which may not be within experimental reach. 

VI. DISCUSSION AND CONCLUSION 

In this paper, we have presented some detailed (analytical 
and numerical) calculations of a model describing a quantum 
dot embedded in two short wires which are connected to reser- 
voirs. We have shown that finite size effects may occur when 
the Kondo temperature T K (for infinite wires) becomes of or- 
der of the level spacing inside the wires. The main conse- 
quence is that the genuine Kondo temperature of this system 
Tk may differ considerably from T K and furthermore be- 
comes sensitive to the local fine structure of the local density 
of states seen by the quantum dot. The LDOS seen by the dot 
can be controlled by appropriate gate voltage . It is also 
worth noting that experimentally it provides an extra parame- 
ter to control the Kondo temperature. This strong difference 
between Tk and T K has a signature in all transport and also 
thermodynamic properties as we have seen. An interesting 
smoking gun to detect these finite size effects is the halving of 
the period (see Figure |4ji in the curve of conductance versus 
wire gate voltage e^ L . 

The symmetric version of the model analyzed in this paper 
has been investigated recently by Cornaglia and Balseiro^i us- 
ing the numerical renormalization group. Their main conclu- 
sions is in agreement with ours (see also Ref. 1 20]). Nev- 
ertheless, on a more quantitative footing, we note that these 
authors found a non-monotonic behavior in some plots of the 
conductance as a function of temperature (particularly for the 
off resonance case). We have not found signatures of such be- 
havior. This might be due to the breakdown of the SBMFT in 
the regime T > T^ R . We also note that a rather large width 
of the peaks in the LDOS was considered by Cornaglia and 
Balseiro (they typically take S ~ T K ) which prevents a dis- 
tinction between the on-resonance case for T K -C A and the 
standard case of infinite wires. 

The analysis developed in the paper relies on some specific 
approximations we want to discuss. 

First, we have assumed that all the conductors, the short 
wires included, are in a ballistic regime, or at least that the 
mean free path Id is the largest length scale by far. If the mean 
free path Id would be for example smaller than the length of 
the short wires, we would then need to compare the Kondo 
length scale ^ with Id and not with d the wire length. 

We have also neglected the Coulomb interactions in the 
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wires. This is one of the most important approximation made 
in our treatment (which also makes the problem tractable). 
One open question is whether (and how) the genuine Kondo 
temperature of the system depends on the Coulomb interac- 
tions in the wires. The Coulomb interactions may also mod- 
ify the transport properties through the whole system. For 
example, Coulomb blockade effects may appear in the con- 
ductance as a function of the wire gate voltage. Depend- 
ing on the local density of states being on or off resonance, 
the conductance through the central quantum dot will be very 
large or intead very small (when finite size effects occur). As- 
suming Coulomb blockade phenomenon are important in the 
wires, such a system can be seen as a double quantum (large) 
dot connected by a weak/strong tunneling junction, a device 
which was studied in Ref. 12911 . When the tunneling junc- 
tion is weak it was found that the period of the peaks doubles 
compared to the strong tunneling limit. Indeed, for the small 
tunneling junction, the excitation spectrum consists of two in- 
dependent spectra of the two dots/wires of size L whereas for 
a strong tunneling junction, the excitation spectrum is bet- 
ter described by the one of a single larger dot of size 2L*& 
Therefore it may be tempting to argue that our prediction of 
the halving of the period in the conductance function of the 
wires gate voltage still occur when Coulomb interactions are 
incorporated. Nevertheless, a complete analysis incorporating 
consistently both the Coulomb interactions and the finite level 
spacing in the wires would be required and goes belong the 
scope of the present paper. We hope to return to this more 
difficult but important problem in the future. 



APPENDIX A: CALCULATION AND ANALYSIS OF THE 
LOCAL DENSITY OF STATES 

In this appendix, we give the main steps of the derivation 
of the expression for the local density of states given by Eq. 
l|9} in our tight binding model. We also study its main charac- 
teristics. In the this appendix, we calculate the left LDOS, the 
right one is obtained by changing L — > R. 

We write the wave function in site j with momentum k as 

tp(j) = Asin(k L j) for - d L < j < -1 

ip(j) = Bsm(kj) + C cos(kj) for j < —dz — 1 

where the wave vectors k and fcj, are related through 

—2t cos k — [i = —2t cos kL + e^ v — fi. 



B, C are related to A by writing the Schroedinger equations at sites — dL and —dL — 1: 



B 



A cos(fcdi) sm(k L (d L + 1)) - 7^ cos(fc(dL + 1)) sin(fc L di,) 



It 



sin k 



(Al) 
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_ A - sin(kd L ) sm(k L (d L + 1)) + 7 2 sm(k(d L + 1)) s'm(k L d L ) 

t> — : — ; , (A4) 

7L S111K 

for 7/, = t WL /t 7^ 0. Moreover, suppose our system is embedded in a huge box of size L » d. Normalizing the wave 
functions to 1 in the box gives the relation \B\ 2 + \C\ 2 — 2/L. Therefore, A 2 reads 

, 9 2 7? sin 2 k 

j^- — . 'J± (A3) 

L sm 2 (k L (d L + 1)) - 2jI cos k sin(fc Z/ ((iL + 1)) sm(fc L d L ) + 7^ sin 2 ^ d^) 

Notice that for Jl = 1 and e w = 0, we recover A 2 = 2/L. The local density of states is then given by pl{c) = 
A 2 sin 2 kL / (2t s\xik(n / L)) with our normalization of the anti-commutation relations. We therefore obtain as a final expres- 
sion the equation (|9j. 

There is another way to derive such expression which is more convenient in order to study the properties of the local density 
of states. We consider the same geometry described in Figure 1. Instead of the quantum dot, we assume a direct tunneling of 
amplitude t' between the sites —1 and 1. We also assume the geometry to be completely symmetric (i.e dL = dn, e w = e w and 
7l = 7ii)' For small t' the conductance can be calculated perturbatively in t': 

G = ^4(t') 2 J de(n PL (e)) 2 (-dn F (e)/de) + 0(t' 4 ) (A4) 
Since the system is non interacting, the conductance can be also computed exactly using the Laudauer formula: 

G =^kJ de\T tot \ 2 (e)(-dn F (e)/de) (A5) 

where T tot (e) is the total transmission amplitude through the system. For a wave number fe, T(k) reads: 

rp ( , x _ T LW T WL T 

JtotW ~ 1 - 2R Q R WL e 2 ^ + (Rwl) 2 (Rq - T 2 )e 2 ^ (AD} 

where Rwl is the reflection coefficient for the left and right weak link in sites — di and da, Tlw is the transmission coefficient 
through the weak link for an electron coming from the lead, Twl is the opposite and finally Ro,To are the reflection and 
transmission coefficients through the central weak link. Since 

, s -2it'smk L „ s t' 2 - 1 

T °W= 1 _ t * f *kr. i B°(k)= 1 _ v2f p kL (AT) 

the term in t' 2 is easily obtained by replacing in the denominator of Eq. (IA6> Ro(k) by —1 and To(fc) by and keeping the first 
term in T (k) in the numerator. Therefore, we can write 



7rpi(fc) = sin k 



Tlv/Twl 



sink L - tlwtwl — — ( A g) 

1 - 2r wi cos a + 



a = ir + 2k L d L + Arg{R WL ) = 2k L (d L + 1) + arctan - ' L 2 y ' - arctan — ' L 2 \ n — (A9) 



where we have defined Rwl = r w ie l9wL , t w i = \Twl\, Uw = \Tlw\ and 

7| sin(fe + fez,) jI sin(fe - k L ) 

— ^ - — - — arctan ^ f 

1 — 72 cos(fc + fez,) 1 — 7£ cos(fc — kL) 

One can check explicitly that the expression in Eq. (IA8> coincides with the one in Eq. l|9}. Fortunately, the equation in JA8> is 
much more compact and readable. One can for example immediately infer from JA9> that the positions of the peaks in the local 
density of states are given by the solutions of a — 2-nn which gives the equation (II 1> . The peak half-width Sl.u can also be 
computed straightforwardly using the expressions 



/ (1 - 7 2 cos(fc - fc £ )) 2 + t| sin (fe - k L ) _ _ 2 lL sin k L 

V (1 - Tl cos(fc + k L )) 2 + H sin (fe + k L ) ^(1 - 7 2 cos(fc + k L )) 2 + 7! sin 2 (fc + k L ) 

The expression for t w i is obtained from ti w by exchanging k with k^. The height of a given peak labelled by n is given by 

/ 7 \ '7 ^IW^WI / A 1 1 \ 

KpL\k n ) = smk L , n j- — (All) 

(*t Twi) 

The width of the peak is obtained for 

6a « 2(d L + l)Sk L = (1 - r wi )/^r^ 
which provides the expression Jl 31 for the width (using Se — 2 sin fci^fci). 
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APPENDIX B: CONDUCTANCE THROUGH A QUANTUM DOT USING PERTURBATION THEORY 



In this appendix, we sketch the derivation of the equation ( I46t . We follow Ref. 12511 and first derive a perturbative expression 
for the current assuming there is a small difference of potential eV — fi^ — fi R between the left and right reservoir. The current 
operator between the left and right side of the dot can be expressed as: 

/(*) = j J dede' r L (e)f R (e') • Sc R , e , + V LR c\ t c R ^ - H.c\ (Bl) 

The current (I) can be for example calculated using 

(I) = (0|S(-oo, 0)/(0)S'(0, -oo)|0) (B2) 

with 



-i J dtH lnt (t) 

S(0, — oo) = Txe "°° 

where Hint contains the Kondo couplings plus the potential scattering terms (see Eq. (I15» . Tk indicates the time ordering 
operator along a Keldysh contour and |0) designs the ground state at time t = — oo. At second order, we simply pick up the 
commutator between the current operator and the interacting Hamiltonian 

= pi? e | J dede'f L {e)f R {e') J (0| Jlrc^O)^ ■ Sc^'(O) + V LR c^ e (<S)c R , e ,{<S), H int (t) |0) J (B3) 

Only the Jl R and Vlr terms give contributions. After performing the Wick contractions and taking the small V limit to extract 
the linear conductance, we recover the second order contributions to the expression in Eq. (I15> . 

The calculation of the third order terms goes along the same line. By expanding both S operators to second order, we obtain 
three sorts of terms: a cross term mixing the first order expansions of S(— oo, 0) and S(0, —oo) (we denote this term l[ 3 ^) and 

(3) 

two other terms coming from the second order expansion of 5(0, — oo) (denoted I\ ) and from the second order expansion of 
S(-oo, 0) (denoted 4 3) ). 

Focusing first on the Kondo terms, the only possible contributions at third order are in J\ r Jll and J\ R J RR . Let us calculate 
for example the terms in J\ R J RR . The mixed term in J RR reads: 



4%= -^JIrJbrJ J j dede'dwIm{p L {e)p R {e')p R {u>) j J ' dhdt 2 {S a (h)S b (0)S c (t 2 ))Tr (yyy 

((c i ,e(ti)4 )£ (0))(c H , e ,(0)4 je ,(t 2 ))(4 ia; (ii)c R)a; (t 2 )) (B4) 

+<4, e , {tx)c R , e , (0)) (4 )6 (0)c iie (*a)) (cr,*, (h)c RiU fa))) } 

Note that no time ordering is necessary here since the crossed terms belong to two different branches of the contour. Using 
(S a (t 1 )s b (0)s c (t 2 )) = |e abc and Tr f^^^J = |e abc (where e abc is the antisymmetric unit tensor), and performing the two 
time integrations, we can easily evaluate this expression. The final result reads: 

R\ 



where n L / R is the occupation number in the left/right reservoir and PP stands for principal part. 

We can perform similar calculations for the two other terms denoted 1^ R and 1$ r paying attention to the different time 
ordering. We give directly the final expressions for these two terms: 
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and 

(3) 2e67r 2 f [,J pp pW(Of n M « n fl A 

+ pp pWffi nf _ nL R _ 1} \ (B7) 
e-e' J 

Adding this three contributions, we find that the terms in p L (e)(p R (e')) 2 cancel out. The terms in J 2 r Jll can be read out 
directly by exchanging L <-* R. For small voltage, we use n L (e) — n R (e) — eV(—dn/de) which provides the linear conductance 
given in Eq. (I15> . 

A priori, one may expect other contributions at third order. The terms in V 2 r Jll and V 2 r Jrr are trivially 0. One can also 
prove by a similar calculation that the terms in J\rVll and J^rVrr involve no logarithm divergencies but instead the integral 
PP J de'p L > R (e')/(e — e'), an integral which is zero in the continuum limit with Do — ► oo. The pure scattering contributions 
in VlrVll and VI r Vrr involve the same integral and therefore do not contain logarithmic divergences as it should be. These 
terms are negligible in the infrared limit since they do not renormalize as oppposed to the terms involving only the Kondo 
couplings. 



APPENDIX C: CONDUCTANCE THROUGH A RESONANT NON INTERACTING QUANTUM DOT CONNECTED TO FINITE 

SIZE WIRES 

In this appendix, we calculate the conductance through the device depicted in Figure 1 when the quantum dot is modeled by 
a non interacting resonant level. It corresponds to U = in (|3j. Generalizing (IA6> to an asymmetric geometry, the transmission 
amplitude through the system 

rplr rplr rplr 

rplr _ 1 LW 1 WL 1 a 

tot\ ) i _ 7-,; nr p 2ik L d L _ f?r tdI p 2ik R d R _|_ ftr rtl I r>l r>r _ rpl r rp r l\ 2ik L d L +k R d R V ; 

where R r LW is reflection coefficient for the left weak link (connecting the left lead to the wire) for a wave coming from the 
right. R l WL is the opposite. T l £ w is the transmission coefficient through the left weak link for an electron coming from the lead 
(from left to right), T^ L is the same for the right weak link, and finally Rq, R 1 ,Tq , Tq 1 are the 4 reflection and transmission 
coefficients through the central weak link. We write 



with 



and 



Tq 1 " = t e lVa+3i ^^ 2 ; R r = r Q e l{lpo+eo) ; R l = r e l(ipo ~ eo) (C2) 



t 2 = = 4(^ g ) 2 (^) 2 sin 2 k L (c3) 

(e + cos/c L + cosfc K - (t^ D ) 2 cosk L - (tfy D ) 2 cos k R ) + ((t^ D ) 2 smk L + (tfy D ) 2 sin k R ) 
2 _ _ (fo + co s k L + cos k R - (t^ D ) 2 cosk L - (^ g ) 2 cos k R ) 2 + ((t^ D ) 2 sinfc L - (t^,) 2 sinfc fl ) 2 

Tq — — 2 2 V 

(e + cosk L +cosk R ~ (t^ D ) 2 cosk L - (^ D ) 2 cos k R ) + ((t^ D ) 2 smk L + {t^ D ) 2 sin k R ) 



(twp) 2 sin k L + ft^ D ) 2 sin k R 
e + cosk L + cosk R - (t^ D ) 2 cosk L - (i^ D ) 2 cosfc fl 

«„ = au-lau ( ^ D ) 2 sink L -(t R D ) 2 sinkR \ 

e + cosk L + cosk R - {tw D Y cosk L - {t^ D ) 2 cos k R J 

(CI) 



W = arc * an ( , — , Z. WZ A I (C5) 



Using 



2i7isinfce^ fc kL "> r cos kL — 7^ cos k + i(sin fe/, — sin A;) 

cos/cl — 7^ cosfc — i(sin/c_L + t£ sinfc) ' LW/ cos fa, — 7£ cos k — i(sin fei + 7^ sin A;) 



WL 



2ijR sin fcfl e 4 ' fc fe -f0 cos fcij — 7^ cos fc + i(sin k R — 7^ sin k) 

cosk R — j R cosk — i(smk R + "f R smk) ' lVL cos k R — j R cos fc — i(sin k R + j R sin k) ' 



- 7fj cos / 



(C9) 
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and the property R l R r - T^T^ 1 = e 2i ^°, the expression O} can be worked out easily. The final result reads 

\rp ,2 ( 2t olLlR sinfc sin k R ) 2 

\J-tot\ c 2 + D 2 ( 



with 



and 



C = cos(fcz,di + k R d R + ipo) (cos(fci + k R ) + JlIr cos 2k — 7^ cos k R cos k — j R cos k^ cos k) 
-r cos(k R d R + 8 Q - k L d L ) (cos(fc L + k R ) + -flj R - 7I cos k R cos k - j R cos k L cos k) 
— sin(fcL<iL + k R d R + <fo) (sin(fci + k R ) — 7^ sin k R cos k — j R sin fc^ cos fc) (CI 1) 

— ?*o sm(k R d R + 9q — kj^dh) (sin(/c£ — k R ) — 7^ sin k R cos fc + 7^ sin fc^ cos fc) 

D = — cos(fc i d i + /c^dfl + (/3 ) sinfc (7^, cos fc L + 7^ cos fcfl — 27^,7^ cos /c) 
r coa(k R d R + - k L d L ) sin fc (7^ cos fc L - 7I cos k R ) 

+ sin(k L d L + k R d R + (p ) sin k (j R cos k L + ~j\ cos k R ) (C12) 
-r sin(k R d R + 9 - kLdL) sin fc (7^ sin kL - 7I sin k R ) 

I 
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